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1.  Introduction 


Shape  representation  is  an  important  aspect  of  image 


processing.  It  is  useful  to  have  a  representation  that  allows 
easy  calculation  of  shape  features  and  shape  differences . 

This  paper  presents  several  algorithms  that  operate  on  regions 


represented  by  quadtrees  [1],  In  particular,  a  linear  time 


method  of  calculating  the  moments  of  a  region  is  presented,  as 


well  as  algorithms  for  the  operations  of  union,  intersection 


It  is  assumed  that  the  image  is  a  square  binary  2  x2 


array  composed  of  black  and  white  points.  The  quadtree  repre 


sen tation  o 


division  into  quadrants.  The  root  node  of  the  tree  represents 


the  whole  image.  If  this  is  not  all  black  or  all  white,  it  is 


split  into  four  quadrants,  represented  by  four  sons  of  the 


part  of  the  image  that  is  all  white  or  all  black  is  subdivided 


to  four  sons.  Eventually,  this  process  terminates,  at  worst 


with  individual  pixels  as  leaves  in  the  tree 


2.  Definitions 

Each  node  in  the  quadtree  holds  six  pieces  of  information. 
The  first  five  are  pointers  to  the  node’s  father  and  to  its 
four  sons. (The  root  node  has  no  father,  and  leaf  nodes  have 
no  sons.)  The  sixth  piece  of  information  describes  the  node 
itself.  It  has  value  BLACK  if  the  region  represented  by  the 
node  is  all  black,  WHITE  if  it  is  all  white,  and  GRAY  other¬ 
wise.  GRAY  nodes  represent  non-terminals  in  the  quadtree, 
while  BLACK  and  WHITE  nodes  are  leaf  nodes. 

The  four  sons  of  a  node  are  labeled  NW  (northwest) ,  NE 
(northeast) ,  SW  (southwest) ,  and  SE  (southeast) .  Notice  that 
if  a  node  is  of  size  2nx2n  (i.e.,  it  corresponds  to  a  region  of 
2"x2n  pixels  in  the  image),  its  sons  will  each  be  of  size 
2n"1x2n~1. 

The  square  array  of  pixels  on  which  the  quadtree  is  defined 
has  the  following  coordinate  system  imposed  on  it. 

The  leftmost,  topmost  pixel  is  at  the  (0,0)  coordinate. 

The  rightmost,  bottommost  pixel  is  at  the  (2n-l,2n-l)  coor¬ 
dinate.  The  coordinates  (x,y)  are  such  that  the  x-value  in- 
| 

creases  from  left  to  right  (west  to  east) ,  and  the  y  coordinate 
increases  form  top  to  bottom  (north  to  south)  in  the  image. 


The  particular  coordinate  system  chosen  is  not  crucial;  this 
one  makes  the  calculations  particularly  easy. 


3 .  The  algorithms 


All  the  algorithms  to  be  presented  here  are  based  on  a 
very  simple  tree-traversal  algorithm.  In  order  to  compute 
moments,  it  is  necessary  to  visit  every  BLACK  leaf  in  the 
tree  and  calculate  some  function  on  the  basis  of  the  distri¬ 
butions  of  the  leaves.  For  the  set  operations,  comparisons 
between  the  structure  of  two  quadtrees  require  parallel  tra¬ 
versal  of  the  trees. 

Where  the  quadtree  representation  scores  over  alternative 
array-based  representations  is  in  its  ability  to  process  large 
blocks  of  the  image  at  once,  when  these  blocks  form  single 
leaves  in  the  quadtree.  Samet  [2)  has  shown  that  the  worst 
case  size  of  a  quadtree  is  4Bn,  where  B  is  the  number  of  BLACK 
nodes  and  the  image  is  of  size  2n  by  2n.  The  expected  size 
should  be  much  smaller,  so  that  substantial  processing  gains 
may  be  achieved. 

The  tree  traversal  algorithm  uses  a  form  of  pcstcraer  L J ] .  Tne 
procedure  . s  to  traverse  the  NW  subtree,  traverse  the  NE  sub¬ 
tree,  traverse  the  SW  subtree,  traverse  the  SE  subtree,  and 
then  visit  the  node.  It  should  be  noted  that  the  particular  tra¬ 
versal  order  chosen  is  not  crucial  for  the  set-theoretic  algorithms. 
Equivalent  algorithms  could  be  devised  for  the  other  orders. 

Al:  Area 

This  algorithm  finds  the  area  of  an  image  represented  by 
a  quadtree,  where  area  is  defined  as  the  total  number  of  black 


pixels  in  the  image. 


Input  to  the  algorithm  is  a  pointer  to  the  root  node  of 
a  quadtree,  and  a  number  n  denoting  the  log  of  the  diameter 
of  the  image  (i.e.,  the  image  is  of  size  2nx2n) . 

integer  procedure  AREA (QUADTREE, N) ; 

/*  find  the  number  of  black  pixels  in  an  image  of  size 
N  Nf 

2  by  2  represented  by  a  quadtree  •/ 

beg  i  n 

node  QUADTREE; 
integer  BLACKAREA; 
level  N; 
quadrant  I ; 

BLACKAREA :-0; 

if  GRAY (QUADTREE)  then 

for  I  in  { NW, NE , SW , SE }  do 
BLACKAREA : -BLACKAREA* 

AREA (SON (QUADTREE, I ) , N - 1 > ; 
else  if_  BLACK  (QUADTREE)  then 

BLACKAREA  -.-BLACKAREA*  2  *  (2*N)  : 
return (BLACKAREA) ; 

2nd ; 


It  can  easily  be  seen  that  this  algorithm  vists  every  black 
leaf  once  and  only  once,  thus  giving  the  correct  area. 

Algorithms  for  finding  higher  moments  are  complicated 
by  the  necessity  of  relating  the  position  of  a  node  in  the 
quadtree  to  the  coordinates  of  a  region  in  the  image.  Each 
quadrant  is  specified  by  the  coordinates  of  its  top  left 


corner,  and  by  its  size.  We  give  an  algorithm  for  finding  the 
centroid  of  an  image.  The  other  moments  are  computed  similarly. 

A2 ;  Centroid 

The  centroid  of  a  binary  image  is  a  point  (x,y)  such  that 

x  is  the  weighted  average  value  of  the  x-coordinates  of  all 

the  black  points  of  the  image  and  y  is  the  weighted  average 

of  the  y-coordinates  of  the  black  points.  In  other  words,  if 

there  are  m  black  points  in  the  image,  (x^ ,y^) , . . . , (xm»ym) » 

“Yi 

the  centroid  is  (x,y)  *  ( - ,  — — )  . 

m  m 

The  centroid  procedure  is  called  with  a  pointer  to  the  root 
node  of  the  quadtree  as  the  QUADTREE  parameter.  N  is  initially 
ecual  to  n,  where  the  input  image  is  size  2n  by  2“;  all  the 
other  parameters  have  value  0. 

procedure  CENTROID (QUADTREE , N , XCOORD , YCOORD , X , Y , MASS , XCENT , YCENT) ; 

N  N 

/*  calculate  the  centroid  of  QUADTREE  for  an  image  of  size  2' x2‘ . 
XCCORD ,  YCOORD  are  the  coordinates  of  the  top 
left  corner  of  the  tree. 

X , Y  are  the  coordinates  of  the  centroid. 

MASS  is  the  number  of  points. 

XCENT  and  YCENT  are  the  centroid  values.  */ 
begin 

node  QUADTREE; 

coord  X , Y , XCOORD , YCOORD , X 1 , Y 1 , XCENT , YCENT ; 
level  N ; 

integer  MASS ,M, J , K . Ml ; 
quadrant  I ; 


X:-0:  Y :  ■  0  ;  MASS:-0; 


if  GRAY (QUADTREE)  then 


for  (I,J,K)  in  { <NW,0,0) ,  (NE,1,0) ,  (SW,0,1)  ,  (SE,1,1)  }  do 
begin 

CENTROID  (SON  (QUADTREE,  I)  ,  N- 1 ,  XCOORD+  J*  2  t  (N-l)  , 
YCOORD+K*  2  t (N-l) , Xl , Y1 ,Ml , XCENT , YCENT)  ; 

X : -X+Xl ; 

Y:-Y+Yl; 

MASS :-MASS+Ml; 
end; 

else  if  BLACK (QUADTREE)  then 
begin 

X:-XCOORD+2t (N-l) ; 

Y:«YCOORD+2f (N-l)  ; 

MAS S  :  “ 2  *  (N+ 1 ) 
end ; 

if  MASS-0  then 
begin 
XCENT :-0; 

YCENT: -0; 
end; 

else  begin 

XCENT : -X/MASS ; 

YCENT : -Y/MASS ; 
end; 


eturn  : 


The  main  difference  between  the  centroid  and  area  algo¬ 
rithms  is  in  the  calculation  of  the  coordinates.  This  is  a 
simple  process,  and  does  not  affect  the  rest  of  the  algorithm. 
It  is  easy  to  see,  once  again,  that  each  BLACK  leaf  in  the 
tree  is  visited  once  and  only  once. 

The  other  moments  can  be  calculated  in  an  analogous  way. 
Set  Operations 

It  is  often  desired  to  compare  two  regions  in  an  image  and 
to  find  what  is  common  to  them.  This  involves  intersecting 
the  two  regions  and  determining  what  points  are  common 
to  both  of  them.  Other  set  operations  that  are  useful  are 
complementation  and  union.  The  algorithms  for  these  operations 
involve,  once  again,  a  tree  traversal.  However,  in  the  union 
and  intersection  algorithms,  it  may  not  be  necessary  to  tra¬ 
verse  the  whole  tree.  The  union  and  intersection  algorithms 
are  special  cases  of  the  superposition  algorithm  of  Hunter  and 
Steiglitt  [ 4 ! .  They  are  included  because  of  their  simplicity 
and  their  similarity  to  the  other  algorithms. 

A3  iComplement 

Constructing  the  complement  of  an  image  involves  changing 
black  pixels  into  white,  and  white  pixels  or  nodes  to  black. 
This  is  a  very  simple  operation  in  a  quadtree,  and  does  not 
change  the  structure  of  the  tree  at  all. 

Input  to  the  procedure  is  a  pointer  to  the  root  of  the 


tree . 


procedure  COMPLEMENT (QUADTREE); 


/*  change  a  quadtree  into  its  complement  */ 
begin 

node  QUADTREE; 
quadrant  I; 

if  GRAY (QUADTREE)  then 

for  I  in  {NW,NE,SW,SE;  do 

COMPLEMENT (SON (QUADTREE, I) ) ; 
else  if_  BLACK  (QUADTREE)  then 
TYPE (QUADTREE) : -WHITE: 
else  /*  a  WHITE  node  */ 

TYPE (QUADTREE) ; -BLACK; 

end; 

A-< ;  Intersection 

This  procedure  finds  the  logical  AND  of  two  binary  images 
represented  by  quadtrees.  The  algorithm  involves  traversing 
the  trees  m  parallel-  When  one  tree  has  a  son  that  is  a  BLACK 
leaf,  while  the  other  has  a  corresponding  son  that  is  not 
3 LACK ,  the  3LACK  leaf  is  replaced  by  the  corresponding  subtree. 
If  one  tree  has  a  leaf  that  is  WHITE,  the  intersection  tree 
will  have  a  corresponding  WHITE  leaf.  Finally,  if  both  trees 
have  GRAY  nodes  in  corresponding  positions,  the  nodes'  sens 
are  examined  recursively,  using  the  same  process. 

Input  to  the  procedure  is  a  pointer  to  the  root  of  each 


quadtree  procedure  INTERSECTION (TREEl , TREE2 ) ; 

/•  returns  the  intersection  of  TREE1  and  TREE2  */ 
begin 

node  TREEl, TREE2, INTERSECT; 
quadrant  I ; 

if  BLACK (TREEl)  0£  WHITE (TREE2)  then 
return (COPY (TREE2 ) ) ; 

else  i_f  BLACK  (TREE2)  or  WHITE  (TREEl)  then 
return (COPY (TREEl) ) ; 

INTERSECT: -CREATENODE () ;  /‘create  a  root  node  */ 
for  I  in  NW,  NE ,  SW,  SE }  do 
begin 

SON (INTERSECT, I! :» INTERSECTION (SON (TREEl ,1) ,SCN (TREE2 ,1) ) 
FATHER ( SON ( INTERSECT ,!)):■ INTERSECT ; 

end : 


return (INTERSECT) ; 


quadtree  procedure  COPY (TREE) ; 

/•  creates  a  tree  structure  identical  to  TREE  */ 
begin 

quadtree  TREE , NEWTREE ; 
quadrant  I ; 

NEWTREE : -CREATENODE ( ) ; 

/*  create  a  node  with  NULL  FATHER,  SON,  and  TYPE  nodes 
TYPE (NEWTREE) : -TYPE (TREE) ; 
for  I  in  {NW, NE.SW.SE,}  do 

lf_  SON  (TREE.  I ) -NULL  then  SON  (NEWTREE  ,  I )  :«NULL; 
else  begin 

SON (NEWTREE, I) : -COPY (SON (TREE , I ) )  ; 

FATHER ( SON (NEWTREE , I ) ) : -NEWTREE ; 
end: 


eturr.  NEWTREE )  ; 


A5 : Union 

The  union  algorithm  is  very  similar  to  the  intersection 
algorithm.  The  trees  are  again  traversed  in  parallel,  and 
decisions  are  made  whenever  either  of  the  traversals  reaches 
a  leaf.  The  decisions  are  the  mirror  images  of  those  in  the 
intersection  algorithm:  when  a  BLACK  leaf  is  encountered,  this 
becomes  the  subtree  at  the  current  position  in  the  union  tree. 
A  WHITE  leaf  in  one  tree  results  in  the  corresponding  subtree 
of  the  other  tree  becoming  the  subtree  at  the  current  position 
of  the  union  tree.  When  there  are  two  GRAY  nodes,  the  proce¬ 
dure  is  called  recursively  for  the  subtrees  rooted  at  these 
nodes . 

Input  to  the  algorithm  is  a  pointer  to  the  root  node  of 


each  tree. 


quadtree  procedure  UNION (TREEl, TREE2) ; 


/*  construct  the  union  of  TREEl  and  TREE2  */ 
begin 

node  TREEl , TREE 2 , UNI ; 
quadrant  I ; 

i_f  BLACK  (TREE2)  0£  WHITE  (TREEl)  then 
return (COPY (TREE2 ) ) 

else  if  WHITE (TREE2)  or  BLACK (TREEl)  then 
return (COPY (TREEl) )  ; 

UNI : "CREATENCDE ( ) ;  /*  create  the  root  node 

for  I  in  {NW,NE,SW,SE}  do 
begin 

SON (UNI, I) : “UNION (SON (TREEl , I) , 

SON (TREE2 , I) ) ; 
FATHER (SON (UNI, I) ) :«UNI; 

end: 

return (UNI )  ; 


Both  the  union  and  intersection  algorithms  can  be  gener¬ 
alized  to  deal  with  an  arbitrary  number  of  trees.  In  all 
cases,  the  intersection  algorithm  takes  time  proportional  to 
the  traversal  time  of  the  smallest  tree,  and  the  union  algo¬ 
rithm  takes  time  proportional  to  that  of  the  second  largest 
tree.  All  the  algorithms  presented  have  times  linear  in  the 
sizes  of  the  trees  on  which  they  operate. 
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»his  paper  describes  algorithms  for  computing  geometric  proper¬ 
ties  of  binary  images  represented  as  quadtrees.  All  the  algo¬ 
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